Exotic 5c-like molecules in QCD sum rules 
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We use the QCD sum rules to study possible iJ^-like molecular states. We consider isoscalar j'' = 0^ and 
J'' = 1^ £)'*'B**' molecular currents. We consider the contributions of condensates up to dimension eight and 
we work at leading order in a.,. We obtain for these states masses around 7 GeV. 
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The study of states with configuration more complex than 
the conventional qq meson and qqq baryon is quite old and, 
despite decades of progress, no exotic hadron has been con- 
clusively identified. Famous examples of possible noncon- 
ventional meson states are the light scalars and the X(3872) 
Jill. While, from a theoretical point of view the most accept- 
able structure for the light scalars is a tetraquark (diquark- 
antidiquark) configuration f2J\, in the case of the X(3872) there 
is an agreement in the community that it might be a DD* 
molecular state. Establishing the structure of these states and 
identifying other possible exotic states represents a remark- 
able progress in hadron physics. 

Besides the X(3872), in the past decade, more and more 
charmonium-like or bottomonium-like states were observed 
in the e+e" collision oI-dII, B meson decays lllla-lsl] and even 
yy fusion processes ll9l4llll. which have stimulated the ex- 
tensive discussion of exotic hadron configurations (for a re- 
view see Refs. Ill2l4l5ll ). An important question that aiTises 
is that if some of these observed states are molecular states, 
then many others should also exist. In a very recent pub- 
lication ilal . a one boson exchange (OBE) model was used 
to investigate hadronic molecules with both open charm and 
open bottom. These new structures were labelled as Bf-like 
molecules, and were categorized into four groups: DS, £)*S*, 
D*!B and DS*, where these symbols represent the group of 
states; D*** - [D'*'",D^*'+,D^*'+] for charmed mesons and 
g(*) = [5W+^ 5(»)0^ 5(*)0] fQj. bottom mesons. A complete 
analysis, based on the approach developed in Refs. Ill7l - t24ll . 
was done in Ref. Illql to study the interaction of these be- 
like molecules. These states were categorized using a hand- 
waving notation, with five-stars, four-stars, etc. A five-star 
state implies that a loosely molecular state probably exists. 
They find five five-star states, all of them isosinglets in the 
light sector, with no strange quarks. 

Here we use the QCD sum rules (QCDSR) 1141125^2711, to 
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check if some of the five-star states found in Ref. llTql are sup- 
ported by a QCDSR calculation. The states we will consider 
are the isosinglets j'' = 0+ DB = (D"B+ + D+B"), 7^ = 1 + 
D*B = (D*"B+ + D*+B"), y^ = 1+ DB* = (D"B*+ + D+B'°) 
and the J^ = 0+ D*B* = (D*°B*+ + D*+B"). The QCDSR 
approach is based on the two-point coiTelation function 



n{q) 



-'! 



d\x e'' 



-'XO|r[;(x)/'"(0)]|0>, 



(1) 



where the current j{x) contains all the information about the 
hadron of interest, like quantum numbers, quarks contents 
and so on. Possible currents for the states described above 
are given in Table HI where we have used a short notation 
for the isoscalars since we are considering the light quarks, 
q - u, ^, degenerate. We use the same techniques developed 
in Refs. 1128143911. 



TABLE I: Currents describing possible Sc-like molecules. 



State 


iW) 


Current 


DB 


0(0+) 


j = (qyscKbysq) 


D'B' 


0(0+) 


j = {qj,c)(by''q) 


D'B 


o(r) 


jfi = i{qy^c)(byf,q) 


DB' 


0(1+) 


j,, = i(qjsc)(hj^q) 



The QCD sum rule is obtained by evaluating the correla- 
tion function in Eq. ([U in two ways; in the OPE side, we 
calculate the correlation function at the quark level in terms 
of quark and gluon fields. We work at leading order in a, in 
the operators, we consider the contributions from condensates 
up to dimension eight. In the phenomenological side, the cor- 
relation function is calculated by inserting intermediate states 
for the hadronic state, H, and parameterizing the coupling of 
these states to the current jftix), in terms of a generic coupling 
parameter A, so that; 



mm) = A, 



for the scalar states and 



{0\j^\H) = A s^. 



(2) 



(3) 



for the axial currents, where e^ is the polarization vector In 
the case of the axial current, we can write the correlation func- 
tion in Eq. ([T]) in terms of two independent Lorentz structures: 



n^M) = -n(r)(^, 
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The two invariant functions, 11 and E', appearing in Eq. (|4|i, 
have respectively the quantum numbers of the spin 1 and 
mesons. Therefore, we choose to work with the Lorentz struc- 
ture gfiv, since it projects out the 1 ^ state. 

The phenomenological side of Eq. ([T]), in the g^y structure 
in the case of the axial currents, can be written as 
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where M„ is the hadron mass and the second term in the RHS 
of Eq. (|5]i denotes the contribution of the continuum of the 
states with the same quantum numbers as the current. In gen- 
eral, in the QCDSR method it is assumed that the continuum 
contribution to the spectral density, p"'"'{s) in Eq. (O, van- 
ishes below a certain continuum threshold sq. Above this 
threshold, it is given by the result obtained in the OPE side. 
Therefore, one uses the ansatz [40| 



p''""{s)=p^''^(s)&{s-so). 



(6) 



The correlation function in the OPE side can be written as 
a dispersion relation: 
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where p'^''^{s) is given by the imaginary part of the correlation 
function: 7rpOra(i) = Im[n''^^(5)]. 

After transferring the continuum contribution to the OPE 
side, and performing a Borel transform, the sum rule can be 
written as 



]2-M 



So 
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where we have introduced the Borel parameter r - l/M^, 
with M being the Borel mass. To extract M„ we take the 
derivative of Eq. (O with respect to Borel parameter r and 
divide the result by Eq. ([8]l, so that: 



M' = 



/ dsse-''^p''P''(s) 
/ dse-'^ p-'P^is) 



(9) 



The expressions forp'''"'(s) for the currents in Table |I] using 
factorization hypothesis, up to dimension-eight condensates, 
are given in appendix lAl 



To extract reliable results from the sum rule, it is necessary 
to establish the Borel window. A valid sum rule exists when 
one can find a Borel window where there are a OPE conver- 
gence, a T-stability and the dominance of the pole contribu- 
tion. The maximum value of t parameter is determined by 
imposing that the contribution of the higher dimension con- 
densate is smaller than 15% of the total contribution. The 
minimum value of r is determined by imposing that the pole 
contribution is equal to the continuum contribution. To guar- 
antee a reliable result extracted from sum rules it is important 
that there is a t stability inside the Borel window. 

The continuum threshold is a physical parameter that 
should be determined from the spectrum of the mesons. Us- 
ing a harmonic-oscillator potential model, it was shown in 
Ref. 1I4III that a constant continuum threshold is a very poor 
approximation. The actual accuracy of the parameters ex- 
tracted from the sum rules improves considerably when using 
a Borel dependent continuum threshold. It also allows to es- 
timate realistic systematic errors IHlll . However, to be able to 
fix the form of the Borel dependent continuum threshold (and 
the values of the parameters in the function) one needs to use 
the experimental value of the mass of the particle [42]. Since 
in our study we do not know the experimental value of the 
masses of the states, it is not possible to fix the Borel depen- 
dent continuum threshold. For this reason, although aware of 
the limitations of the values we are going to extract from the 
sum rule, to have a first estimate for the values of the masses 
of the states, we are going to use a constant continuum thresh- 
old. In many cases, a good approximation for the value of 
the continuum threshold is the value of the mass of the first 
excited state squared. In some known cases, like the p and 
Jjif/, the first excited state has a mass approximately 0.5 GeV 
above the ground state mass. Since here we do not know the 
spectrum for the hadrons studied, we will fix the continuum 
threshold range starting with the smaller value which provides 
a valid Borel window. The optimal choice for .s^o will be taken 
when there is ar- stability inside the Borel window. 

For a consistent comparison with the results obtained for 
the other molecular states using the QCDSR approach, we 
have considered here the same values used for the quark 
masses and condensates as in Refs. ll29l - (35ll43ll . listed in Ta- 
ble ini For the heavy quark masses, we could use the range 
spanned by the running MS mass rngiMQ) and the on-shell 
mass from QCD (spectral) sum rules compiled in ll27ll and 
more recently obtained in Ref. ll44ll . However, we do not ob- 
tain a valid borel window with the usual on-shell mass for b 
quark, m/, = 4.70 GeV. For this reason, we have considered 
as the maximum value for b quark mass mt, = 4.60 GeV, as 
indicated in Table HIl For the (G^) condensate, we have used 
the new numerical value estimated in Ref. ll44ll . To take into 
account the violation of the factorization hypothesis we intro- 
duced in Tablellllthe parameter p. 

Let us consider first the molecular current for the DB (0^) 
state. In Fig. [1] a), we show the relative contribution of the 
terms in the OPE side of the sum rule, for -^^o = 7.20 GeV. 
From this figure we see that the contribution of the dimension- 
8 condensate is smaller than 15% of the total contribution 
for values of t < 0.27 GeV"^, which indicates a good OPE 



TABLE II: QCD input parameters. 



Parameters 


Values 


mt 


(4.24 - 4.60) GeV 


yn. 


(1.23- 1.47) GeV 


im) 


-(0.23 + 0.03)^^ GeV^^ 


{8]G') 


(0.88 ± 0.25) GeV* 


ml = {qGq)l{qq) 


(0.8±0.I)GeV2 


islG') 


(0.58 ±0.18) GeV* 


p = (qqqq)l{qqf 


(0.5 - 2.0) 



a) 




0.15 



0.20 



0.25 
T (GeV^^) 



0.30 



0.35 



convergence. From Fig. [Db), we also see that the pole con- 
tribution is bigger than the continuum contribution only for 
T > 0.22 GeV"^. Therefore, we fix the Borel Window as: 
(0.22 < T < 0.27) GeV"^. The results for the mass are shown 
in Fig. [T] c), as a function of r, for different values of ^o- 
As we can see from Fig. [T]c), the Borel window (indicated 
through the parenthesis) gets smaller as the value of yfsQ de- 
creases. So, we can only work with values for yfsQ bigger than 
7.00 GeV, otherwise we do not obtain a valid Borel window 
for this sum rule. We also observe that the optimal choice 
for the continuum threshold is yfsQ - 7.20 GeV, because it 
provides the best r-stability inside of the Borel window, in- 
cluding the existence of a minimum point for the value of the 
mass. 

Therefore, varying the value of the continuum threshold in 
the range yfsQ — (7.00 - 7.30) GeV, and the others parameters 
as indicated in Table HH we get: 



b) 



m; 



<8> 
DB 



(6.77 + 0.11) GeV. 



(10) 



The quoted uncertainty is the OPE uncertainty. The most im- 
portant source of uncertainty is the values of the heavy quark 
masses. As discussed in Ref. I42[| . there is another kind of 
uncertainty, called systematic uncertainty, related to the in- 
trinsic limited accuracy of the method. The systematic un- 
certainty of the physical quantity extracted from the QCDSR 
represents, perhaps, the most subtle point in the application 
of the method. Without an estimate of the systematic uncer- 
tainty, the numerical value of the physical quantity one reads 
off from the Borel window might differ significantly from its 
true value. In Ref. ll42ll it was shown that the use of the Borel 
dependent continuum threshold allows to estimate the system- 
atic uncertainty. Inparticular, for the case of the D and Ds 
mesons studied in 14211 . the systematic uncertainty turns out to 
be of the same order of the OPE uncertainty. Since here we 
do not have how to estimate the Borel dependent continuum 
threshold, in an attempt to obtain some information about the 
systematic uncertainty, we will repeat the analysis considering 
only terms up to dimension 6 in the OPE. These new results 
are shown in the Fig. |2] 

As one can see in Fig. |2]a), when we remove the dimen- 
sion 8 condensates contribution we lose the OPE convergence, 
since the most important contributions to the OPE come from 
{qq) and p{qq}^ contributions. Thus to be able to extract some 
results from this analysis we determine the maximum value of 



c) 



1.0 



S 0.8 -■■■•■•... 




0.0 



0.15 



0.20 



0.25 
T (GeV^^) 



0.30 



0.35 




FIG. 1: DB (0^) moleculai' state up to dimension 8 contribution for nic = 
1.23 GeV and mt, = 4.24 GeV. a) OPE convergence in the region 0.14 < 
T < 0.36 GeV"^ for ^/so = 7.20 GeV. We plot the relative contributions 
stalling with the perturbative contribution and each other line represents the 
relative contribution after adding of one extra condensate in the expansion: + 
(ijq), + (G*>, + iqGq), + (qq)' and + (qq)(qGq). b) The pole and continuum 
contributions for ^fso = 7.20 GeV. c) The mass as a function of the sum mle 
parameter t, for different values of -\/s(). For each line, the region bounded 
by parenthesis indicates a valid Borel window. 



T parameter imposing that the contribution of the dimension 
6 condensate is smaller than 25% of the total contribution, 
otherwise we do not have a valid Borel window for this sum 
rule. The minimum value of t is not changed since the pole 
dominance behavior remains the same. Finally, we obtain the 
results shown in the Fig. |2]c), from where we get: 



M 



<6> 



= (6.63 ±0.09) GeV. 



(11) 



Note that the value in Eq.(fTTI) dififers at maximum only 
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FIG. 2: DB (0*) moleculai' state up to dimension 6 contribution for m^ = 
1.23 GeV and mj, =4.24GeV. a) OPE convergence in the region 0.19 < r < 
0.27 GeV"" for ^fs^ = 7.20 GeV. We plot the relative contributions starting 
with the peiturbative contiibution and each other line represents the relative 
contribution after adding of one extra condensate in the expansion: + {qq), 
+ (G^), + (qGq) and + (qq)~. b) The pole and continuum contributions for 
■\/so = 7.20 GeV. c) The mass as a function of the sum rule parameter r, 
for different values of y^o. For each hne, the region bounded by parenthesis 
indicates a valid Borel window. 



~5.0% to that in Eq. fTOb . Besides, the inclusion of dimension 
8 condensate provides a better OPE convergence, r-stabihty 
and an improved Borel window. Therefore, even being aware 
that this is only part of the dimension-8 contribution, here we 
consider it as a form to estimate the systematic uncertainty. 
One should note that a complete evaluation of the dimension-8 
contributions require more involved analysis including a non- 
trivial choice of the factorization assumption basis fl45il . Then, 



the final value for the DB molecular state is given by: 
Mdb = (6.75 ±0.14) GeV. 



(12) 



The mass in Eq. ([T2]i is -400 MeV below the DB threshold 
indicating that such molecular state would be tightly bound. 
This result, for the binding energy, is very different than the 
obtained in Ref. iflql for the DB (0^) molecular state. The 
authors of Ref. Il&\ found that the DB (0+) molecular state 
is loosely bound with a binding energy smaller than 14 MeV. 
However, it is very important to notice that since the molecular 
currents given in Table |I] are local, they do not represent ex- 
tended objects, with two mesons separated in space, but rather 
a very compact object with two singlet quark-antiquark pairs. 
Therefore, the result obtained here may suggest that, although 
a loosely bound DB (0^) molecular state can exist, it may not 
be the ground state for a four-quark exotic state with the same 
quantum numbers and quark content. 

Having the hadron mass, we can also evaluate the coupling 
parameter. A, defined in the Eq.©. We get: 



Adb = (0.029 + 0.008) GeV^ 



(13) 



The parameter A gives a measure of the strength of the cou- 
pling between the current and the state. The result in Eq. (|T3T l 
has the same order of magnitude as the coupling obtained for 
the X(3872) ||29I1 . for example. This indicates that such state 
could be very well represented by the respective current in 
Table I] 

We can extend the same analysis to study the others molec- 
ular states presented in Table|I] For all of them we get a similar 
OPE convergence in a region where the pole contribution is 
bigger than the continuum contribution. We obtain the results 
shown in the Fig. [3] 

In Fig. [3] a), we show the ground state mass, for the D*B*, 
0^ molecular current, as a function of t. For -sfso = 7.80 GeV, 
we can fix the Borel window as: (0.18 < t < 0.21) GeV ". 
/,From this figure we again see that there is a very good r- 
stability in the determined Borel window. 

Varying the value of the continuum threshold in the range 
-\/sq - (7.60 - 7.90) GeV, the others parameters as indicated 
in Table HI] and also estimating the uncertainty by neglecting 
the dimension-8 contribution we get: 



Md'B' = (7.27 + 0.12) GeV, 
Ad-b' = (0.115 + 0.021) GeV^ 



(14) 
(15) 



The obtained mass indicates a binding energy of the order of 
~ 50 MeV below the D*B* threshold. Considering the un- 
certainties, it is even possible that this state is not bound. In 
this case, our central result is in a good agreement with the 
result obtained for the D*B*, 0^, molecular state obtained in 
Ref. 11611 . However, since we do not have a trustable estimate 
for the systematic error, as discussed above, any conclusion 
about the possible existence of this state would be premature. 
We now consider the D*B (1^) molecular current. In Fig.|3] 
b), we show the ground state mass, as a function of t. For 
■\fso = 7.70 GeV, we can fix the Borel window as: (0.18 < 
T<0.22)GeV"^ 
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FIG. 3: The mass as a function of the sum rule parameter r, for m^ = 
1.23 GeV and mi, = 4.24 GeV, considering different values for -\/sq: a) for 
D'B', 0* molecular current; b) for D'B, 1+ molecular current; c) for DB', 1* 
molecular current. For each line, the region bounded by parenthesis indicates 
a valid Borel window. 



Varying the value of the continuum threshold in the range 
yfso = (7.50 - 7.80) GeV, the others parameters as indicated 
in Table |II] and also estimating the uncertainty by neglecting 
the dimension-8 contribution we get: 



Md'b = (7.16 ±0.12) GeV, 
Ad-b = (0.058 + 0.013) GeV^ 



(16) 
(17) 



The obtained mass indicates a central binding energy for the 
D*B, 1+ state of the order of ~ 130MeV. Considering the 
uncertainty, this result might be compatible with the one ob- 
tained by the authors in Ref. [|16[1 , or can even be unbound, 
as the D*B* state. Therefore, also in this case, any conclusion 



about the possible existence of this state would be premature. 
Finally, we study the molecular current for DB*, 1^ state. 
As one can see from Fig. [3]c), we have a very good r-stability 
inside of the Borel window: (0.21 < t < 0.25) GeV ^ for 
-\fsQ - 7.30 GeV. Doing the same procedure to estimate the 
uncertainties in the range yfso - (7.10 - 7.40) GeV we get: 

Mdb- = (6.85 ±0.15) GeV, (18) 

Adb' = (0.036 ± 0.01 1) GeV^ (19) 

which indicates a binding energy of the order ~ 330 Me V, 
much bigger than that obtained in Ref. 1II6II . 

We can compare our results with the ones presented in 
Ref. ll46ll . First of all we would like to point out that we have 
found some disagreements in the spectral densities expres- 
sions for the DB, D*B and DB* molecular currents. In par- 
ticular, we have found some missing terms in the (G^), {qGq) 
and {G^) contributions, due to some diagrams that have been 
neglected in their calculations. We have found that the {qGq) 
contribution plays an important role to the final result, and this 
can explain why the mass values found in Ref.|!46«l differ from 
ours. Another important point, in which our calculations dif- 
fer, is the fact that the Borel window (0. 10 < r < 0. 14) GeV"^, 
considered by the authors in Ref. j^ , does not have pole dom- 
inance, as can be seen in Figs. [T]c),[3]a),[3]b) and[3]c). The 
only result for the mass, which is in agreement with Ref. ll46ll , 
is the one for the D*B* molecular current. For this current we 
found disagreements only for the (G^) contribution. Since the 
{G^) contribution is very small, as compared to the others, the 
differences found could not modify the final result. 

In conclusion, we have studied the mass of the exotic Be- 
like molecular states using QCD sum rules. We find that for 
the molecular currents D*B* (7^ = 0+) and D*B{J'' = 1+), the 
QCDSR central results lead approximately to the same pre- 
dictions made by the authors in Ref. 11 1611 . for the respective 
molecular states in a OBE model. However, since our un- 
certainties are underestimated due to our crude model for the 
continuum threshold, any conclusion about the possible exis- 
tence of these states would be premature. 

In the case of the DB (J^ = 0+) and DB* (J^ = 1+) molecu- 
lar currents, from the QCD sum rule point of view, the masses 
of the corresponding states are smaller than the masses ob- 
tained for the respective molecular states studied in Ref. III6I1 . 
We interpret this result as an indication of the possible exis- 
tence of four-quark states, with the same quark content and 
quantum numbers as the DB (j'' - 0+) and DB* (j'' - 1+) 
molecular states, but with smaller masses. 
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Appendix A: Spectral Densities 



For the D*B, 1^ molecular current we get: 



The spectral densities expressions for the molecular cur- 
rents given in Table |I] were calculated up to dimension-6 con- 
densates, at leading order in a^. To keep the heavy quark 
mass finite, we use the momentum-space expression for the 
heavy quark propagator We calculate the light quark part of 
the correlation function in the coordinate-space, and we use 
the Schwinger parameters to evaluate the heavy quark part of 
the correlator To evaluate the d'^x integration in Eq. ([TJ, we 
use again the Schwinger parameters, after a Wick rotation. Fi- 
nally we get integrals in the Schwinger parameters. The result 
of these integrals are given in terms of logarithmic functions, 
from where we extract the spectral densities and the limits of 
the integration. The same technique can be used to evalu- 
ate the condensate contributions. To evaluate the systematic 
uncertainty we also include a part of the dimension-8 con- 
tribution, related with the mixed-condensate times the quark 
condensate. In Ref. pSj it was shown that the contribution of 
this condensate is much bigger than other dimension-8 con- 
densates, related with the gluon condensate. 

For the DB, 0^ molecular current we get: 
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2/2 / 2 2m 

a - m^ - a(m^ - m^) r 



^I'lijov l^fl' 



- 2j^j'-^(J3'm,+a'mh)F{a,/3) 






1^ DR ^ ^ 



a,m„ A 

mi,mcp{qqf 



^|,,2„, , .,2., 



^hc V, 



213 



167r2 

a:,„„r l-a 



^TT^ J ff3j ^3^ 
3 , o3 



2<^ml+aWf,) 



+ {a' +p')F{a,l3) 



1 



P^i^) 



mcmhp{qq){qGq) f da 



2^Tf- 



C da 

Ja(l-a) 



1 — a + ff" 



-\m^ - a(m^ - mf^) )t 



61 s- 



a(l —a) 



For the DB*, 1^ molecular current we get: 



pP"'(s) = 



Pt'^^) 



^i/S/l^'"""^^^'^"^^^"'^"'^^' 



212; 

a,m„ A 

a,,,,, V 1 -a 

3{qq) fda (^d/3 



a,„^n- 1 -a 

m) fda^ fa 

'n^ J a'- jj 

a,„i„ P,„m 



f3{a + P) nic + a nil, 



F{a,(3f, 



<c^>, , <g'G') fda fd/3 
^'DB^ 2i2;r6 J a^ J 0^ ^ 



ami,, P,„ 



al3\3li{a + P) 



a(2-a- P) \F{a,P) + {p^m] + a^ml) 



X {l-a-p)(_l +a+P) 






XqGq) 



28 



ji^ 



J M 



da 



(l-a) 



Oinax l-a 



a,,. 



m.-aOn.-^,))Hia) p<t>(.) = |g> J| J|(l - . -^) 



Mr I — ]d/3{3a + 2p)F(a,p) 



(Xmax i -a 

ami,, P„,i„ 

^™* 167r2 

iC^) Cda rd/3 






G') rda rdp_ 



(l-a-/3)(l +«+/?) 



X 2{0^ml+a^nl) + (a^ +p^)F{a,j3) 



<8> ^ m,mbp{qq){qGq) r da J m--a{m^^-ml) 
^™' " 2^772 jaCl-a) \ a{\-a) 



X 



(1 - af' - m^ - Qf(m^ - mi)\T 



For the D*B*, 0^ molecular current we get: 



<*-/»/<, /?», 



a,„cn 1 -a 



pS!w = -f^/$/|(/i-»,..-».)f(».»=. 






,<G^> 



p' '(s) = 






o3™2 , „,3™2x 



X 08^m^ + a'ml)F{a,fi), 






"»«» A", 



+ (Q'^+yS^^)F(a,j6) 



2(jft7tl + a\nl) 



mcmbp{qq){qGq) f da 



m^-a(m^—nr^)\ 



<8> / N _ "ici'ibhJ\HH/\H^q) C da I ^ "V Liv'"e '"fo/ 



8;r2 



a;(l - a) + I m^ - Q'(mp - m^) )t 



(^(l - a) 



In all these expressions we have used the following defini- 
tions: 



H{a) 
F(a,l3) 

he 



OT^a + m^{\ - a) — a{\ — a)s, 



nif^a + m^./3 - a/3s. 



l+iml- ml)/s, 






and the integration limits are given by: 



(Al) 
(A2) 

(A3) 
(A4) 



A. 



an. 



a,r. 



a m. 



as — m, 

~2 



2' 



(1-v), 
(1+v) . 



(A5) 
(A6) 
(A7) 
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